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ABSTRACT 

Properties of X-ray luminosities in low-mass X-ray binaries (LMXBs) mainly depend on 
donors. We have carried out a detailed study of donors in persistent neutron-star LMXBs 
(PLMXBs) by means of a population synthesis code. PLMXBs with different donors have differ- 
ent formation channels. Our numerical simulations show that more than 90% of PLMXBs have 
main sequence (MS) donors, and PLMXBs with red giant (RG) donors via stellar wind (Wind) 
are negligible. In our model, most of neutron stars (NSs) in PLMXBs with hydrogen-rich donors 
form via core-collapse supernovae, while more than 90% of NSs in PLMXBs with naked helium 
star (He) donors or white dwarf (WD) donors form via an evolution-induced collapse via helium 
star (1.4 < Muc/Mq < 2.5) or an accretion-induced collapses for an accreting ONeMg WD. 

PLMXBs with different donors have different properties. In PLMXBs with MS donors, the 
orbital periods are between ~ 1 hour and 100 hours, and the mass transfer is driven by donor 
evolution or magnetic braking. Our population synthesis code shows that their X-ray luminosities 
mainly are around ~ 10'^® erg s~^. Similarly, in PLMXBs with RG donors via Roche lobe 
overflow(Roche), the mass transfer is driven by donor evolution, but orbital periods are between 
~ 10 hours and 1000 hours. Their X-ray luminosities are around ~ 10'^^ erg s~^. The 2 known 
LMXBs ( Cyg X-2 and GX 13-f 1) can belong to PLMXBs with RG (Roche). PLMXBs with 
RG (Wind) donors have the longest orbital periods and low X-ray luminosities (^ 10^"^ erg s^^). 
Their contributions to X-ray luminosities can be negligible. In PLMXBs with He donors, the 
orbital periods are shorter than 80 minutes, and the mass transfer is mainly driven by magnetic 
braking. Results of our numerical simulations predict that PLMXBs with X-ray luminosities 
around ~ 10'^* erg s~^ mainly come from binaries with He donors. In PLMXBs with WD donors, 
the orbital periods are shorter than 1 hour, and the mass transfer is mainly driven by gravitational 
radiation. According to results of our population synthesis code, their X-ray luminosities are 
between ^ 6 x 10'^^ — 10'^^ erg s~^, and most of LMXBs with WD donors are transient. 

Subject headings: h 



inaries: close — pulsar: general — stars: neutron- 
X-ray: stars 

1. Introduction 

Low-mass X-ray binaries (LMXBs) were discov- 
ered nearly 50 yr ago, and there are n ow ^ 200 
known in the Galaxy (|Liu et all [20071 ). LMXB 



^National Astronomical Observatories / Xinjiang Ob- 
servatory, the Chinese Academy of Sciences, Urumqi, 
830011, China 

^School of Physical Science and Technology, Xinjiang 
University, Urumqi, 830046, China. 



is a mass-transferring binary system with a com- 
pact object accretor(a black hole or a neutron star 
(NS)), and a low-mass (< IMq) donor. The X- 
ray luminosity function is an important character- 
istic of the LMXBs, and has become a key tool for 
studying LMXBs. 

Using results of Chandra observations of old 
stellar systems in 11 nearby galaxies of various 
morphologic al types a n d the census of LMXBs in 
the Galaxy, iGilfanovl ()2004l) suggested that the 
total number of LMXBs and their combined lu- 
minosity are proportional to the stellar mass of 
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the host galaxy. IPostnov fc Kuranov ( 2005 ) sug- 
gested that the flattening of the LMXBs luminos- 
ity function at lower than 2 x 10'^^ erg s^^ might 
correspond to the transition from the magnetic 
stellar wind brakin g to the gravitational wave 
braking mechanism. iRevnivtsev et al.l (|201l[ ) sug- 
gested that LMXBs with X-ray luminosities be- 
low 2 X 10'^'' have unevolved secondary compan- 
ions while systems with higher X-ray luminosity 
predominantly harbor giant donors. 

In theoretical work, people usually assume that 
the X-ray luminosity is directly proportional to the 
mass-accretion rate of compact stars in LMXBs. 
The mass-accretion rate depends on the orbital 
period, donor's evolution, and angular momen- 
tum loss. Donor is usually a main sequence 
(MS) or a white dwarf (WD), and it may also 
be a red giant (RG) or a naked helium star (He). 
The properties of LMXBs are closely related to 
their donors. If donors are poor-hydrogen stars 
(WDs or Hcs), LMXBs usually are ultra-compact 
X-ray binaries wh ose orbital periods are shorter 

Hai). If donors 



2. Models 



than 80 minutes ([Nelson et al 
are RGs 



LMXBs are called as symbiotic X-ray 
Up to now, there 



binaries ( Masetti et al.l 120061 ) 



are about a dozen of ultra-compact X-ray binaries 
and about 10 symbiotic X-ray binaries. 

Different types of stars have different properties 
and evolutions. Therefore, the donors in LMXBs 
determine orbital period and angular momentum 
loss, which directly affects the mass-transfer rate 
from the donors to NSs. Obviously, the donors 
really determine the luminosity of LMXBs. Si- 
multaneously, the observational luminosity func- 
tion of LMXBs provides important constrains for 
our simulating donors' evolution. 

On observations, LMXBs are divided into tran- 
sient and persistent sources. It is difficult to esti- 
mate the X-ray luminosities of transient LMXBs 
during quiescent state and outburst state. In this 
paper, we focus on persistent LMXBs (PLMXBs) 
with accreting NS and different donors, and inves- 
tigate their properties and contributions to total 
X-ray luminosity. In § 2 we present our assump- 
tions and describe some details of the modeling 
algorithm. In § 3 we discuss the main results. In 
§ 4 the main conclusions are given. 



For the simulation of binary evol ution, we use 
rapi d binary star evol ution code BSE (Hurle y et al 



rapi i 
[200^ 



2002,) with updates bv lKiel fc Hurley (20061 BSE 
code calculates the orbital changes of binary sys- 
tems via mass variations, gravitational radia- 
tion and ma g netic braking. Details are in §2 of 
Hurley et all (|2002f) . 



Forming channels of LMXBs are important 
questions in X-ray astronomy. T here are many lit 
eratures to investigate them (e.g . 



1991 



Lin et al 



II2OO2 



■Bhattacharva fc van den Heuvell 
Pfahl et al.ll2'003[ 



Podsiadlowski et al, 

2OIOI) . The difficulty in a theoretical 



investigation involves two problems: (i)keeping the 
binary bound when the massive progenitor of NS 
explodes in a supernova event;(ii)a common en- 
velope (CE) phase. They have great effects on 
LMXBs. The following two subsections give de- 
scriptions. 

2.1. Common Envelope Evolution 

In binary system, due to orbital angular mo- 
mentum loss or stellar expansion, a star can over- 
flow its Roche lobe. If the mass ratio of the com- 
ponents {q = Mdonor/-/V^accrctor) at the OUSCt of 

Roche lobe overflow is larger than a certain critical 
value Qc , the mass transfer is dynamically unstable 
and results in the formation of a CE. The issue of 
the criterion for dynamically unstable Roche lobe 
overflo w, Qc, is stUl open . Based on the polytropic 
models, IWebbinkI (|l988[ ) gave Qc for red giants by 



gc = 0.362 



1 



3 X (1 - Afc/Mdonor) ' 



(1) 



where Mc is the core mass of donor. However, 
this Qc is obtained u n der conserv ative Roche lobe 
overflow. iHan et al.l (|200ll l2002f ) showed that Qc 
depends heavily on the assumed mass-transfer ef- 
ficiency. They found that almost linearly in- 
creases with the amount of the mas s and mo- 
mentu m lost during mass transfer. In iHan et al 



()2002l ) the critical mass ratio qc is between 1.1 and 
1.3. In this work, we take qc —Eq.^, qc = 1.2 and 
qc — 2.0 in different simulations, respectively. 

Although many efforts have been devoted 
to understanding the evolution of CE (e.g.. 



Ricker fc Taaml2008t[Ge et al.ll2010l : rDelove fc Taam 
20101 ). the knowledge about it is still poor. It 

is generally assumed that the orbital energy of 
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the binary is used to expel the envelope of the 
donor with an efheiency ace , which is called as 



a-algorithm. iNelemans et al.l ^2000) suggested to 



describe the CE evolution by an algorithm based 
on the equation for the system orbital angular 
momentum balance wh ich implicitly a ssumes the 
conservation of energ v (IWebbinkl | l984ll. which is 



called as 7-algorithm. Following iLii et alJ (|2006l ). 



for CE evolution in different simulations we use 
ctceAcc = 1-0 in a-algorithm and 7 = 1.5 in 7- 
algorithm, respectively. Here Ace is a structure pa- 
rameter which depends on the evolutionary stage 
of the donor. 

2.2. Formation Channels of Neutron Stars 
and Kick Velocity 

In X-ray binari es, NSs can be formed via 
three channels(e.g., Ivanova et al. 2008HKiel et al. 
20081) : (i) Core-collapse supernovae (CCSN) for 



a star with main-sequence mass M/Mq > 11; 
(ii)Evolution induced collapse (EIC) of a he- 
lium star with a mass between 1.4 and 2.5Mq 
in which the collapse is t riggered by electron cap- 



ture on 20Ne and ^^Mg (jMivaii et al.lll980t ): ( 



Accretion-induced collapses (AIC) for an accret- 
ing ONeMg WD whose mass reaches the Chan- 
drasekhar limit. Response of accreting ONeMg 
WD is treated in the sai ne way as the evolution 
of CO WD (see details in lLii et al.|[2009l ). 

Nascent NS receives additional velocity ( "kick" ) 
due to some still unclear process that disrupts 
spherical symmetry during the collapse or later 
Dichotomous nature of kick s which was suggested 
quite early bv iKatz (1975). Observationally, the 
kick is not well constrained due to numerous selec- 
tion effects. Currently, high kicks 100km s~^) 
are associated with NS originating from CCSN, 
while low kicks (~ 10km s~^) with NS born in 



EIC and AIC (|Pfahl et al.ll2002D . 



We apply to core-collapse NS Maxwellian dis- 
tribution of kick velocity Uk 



2 V 



3^ 



(2) 



CTk = 190 and 400 km s^^ for CCSN, while ajj = 
20 and 10 km s^^ for EIC and AIC in different 
simulations. 



2.3. X-ray Luminosity 

The X-ray luminosity of the accreting NS can 
be approximated by 



iboi - vM^sc^ ^ 5.7 X lO^Serg s-\^){-^ 



Mns 



lO-ii'Mgyr-i ' 

(3) 

where 77 ~ 0.1 is the efficiency of accretion onto 
the NS and Mns is the mass-accretion rate of the 
NS. Super-Eddington accretion rates may be im- 
portant in the formation of low-mass X-ray bina- 
ries a nd millisecond pulsars! Webbink fc Kalogeral 
19971 ). We assume that Mns = niin(MNs, '7Edd x MEdd) 
where MEdd is Eddington limit given by 

A/Edd - 2.08 X 10-3(1 + X)-ii?NsM0yr-i. (4) 

Here, X is the hydrogen mass fraction. 77Edd is 
the factor to al low super-Eddin g ton luminosities , 
taken to be 5 (|Begelmanl |2003 IZuo fc Lil I2OIII ). 
To transform the bolometric luminosity into the 
X-ray luminosity, a bolometric correction factor 
77boi is introduced by L x = ??boiiboi- Following 
Belczvnski et al.1 (|2008l ). we take ryboi = 0.55. 

Roche overflow-fed systems are subject to a 
thermal disk instability and may appear either 
as persistent or transient X-ray sources depend- 
ing on the mass transfer rate. A system becomes 
a transient X-ray source when the mass-transfer 
rate falls below a certain critical value, Afcrit. For 
hydrogen-rich disks ( The donors are MS s or RGs 
), we use the work of Ivan Paradijsl ( 1996f l. Apply- 
ing to Eq.(l3]), Afcrit for hydrogen- rich disk is given 

by 

Merit ^ 1-8 X 10^^ P^' g/s, (5) 

where Porb is orbital period in hours. For 
disks with heavier elements, we use the work of 
Menou et all (|2002l ): 



5.9 X lOi^Af^s' -R: 



-0.87 7j2. 62^,0. 44 



0.1 



?/s, 



He rich 



Merit = 



1.2 X 1016 Afj^O-^4^2.21^0.42 g/g^ crich 

5.0 X 10i6AfNs° ''i?^ °^aO:f g/s, O rich 



(6) 

where Rd is a maximum disk radius (2/3 of accre- 
tor Roche lobe radius) in 10^° cm, ao.i = a/0.1 in 
which a = 0.1 is a viscosity parameter. 

If A/ns > Afcrit or wind-fed accretion, the sys- 
tem is a PLMXB whose X-ray luminosity is deter- 
mined by Eq. ([3]). If AfNs < Afcrit in Roche over- 
flow accretion, the system is a transient source. In 
this work we focus on PLMXBs. 
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3. Results 

We use Monte Carlo method to simulate the 
initial binaries. For initial mass function, mass- 
ratios, and separations of components in binary 
sy stems, we adopt the distributions used by us 
in iLii et al.l (|2006l 120081 ). We assume that all bi- 
naries have initially circular orbits. After a su- 
pernova, new parameters of the orb it are derived 
using standard formulae, (e 



Hameurv et~allll993l : iHurlev et al.ll2010l ). In this 



Hurlev et al 



2002[ ). Table [T] lists all cases considered in the 
present work. Our model is normalized to forma- 
tion of one binary wit h Mi > O.SMq per year 



([Yungelson et al.l 119931 ). We use 1 x lO" binary 



systems in our Monte-Carlo simulations. 

In this work, a binary is considered as PLMXB 
if it satisfies the following conditions: 

(i) Binary includes an NS and its companion's mass 
is lower than 6 Mq; 

(ii) Binary orbital period and the mass-accretion 
rate of NS satisfy Eqs.© and dH). 

Here, we call both low- and intermediate-mass X- 
ray binaries as LMXBs. 



3.1. 



Birthrates and Numbers of PLMXBs' 
Populations 



In 



our simulations, there are 29000 (case 
2) — 110000 (case 3) PLMXBs in the Galaxy, 
and their birthrates are ~ 3.4 — 7.2 xl0~'* yr~^. 
However, the nu mber of ah obs er ved LMXBs 
i s less than 200 (|Liu et al.l bood ). IPfahl et al. 
(|2003 ^ investigated LMXBs via CCSN, and ob- 
tained birthrates for LMXBs of lO'^— lO"'' yr"! 
and 400—70000 LMXBs, a factor 10—1000 times 
higher than observed number. Many authors 
suspected that the mismatch between the ob- 
served number and theoretically predicted num- 
ber could be related to irradiation effects (et al. 



work, we do not consider LMXBs in low states 
driven by irradiation-driven limit cycles. How- 
ever, we still encounter the known problem of 
overproduction of LMXBs. 

As Table [2] shows, NSs with different kinds of 
donors in PLMXBs have different formation chan- 
nels. More than 90% of PLMXBs have undergone 
CCSN, especially, for PLMXBs with hydrogen- 
rich donors. However, more than 90% of PLMXBs 
with He donors or WD donors have undergone 
AIC and EIC. In cases 1 and 2, parameter cTk is 
increased from 190 to 400 km s~^. The larger 
(Tk is, the more difficultly binary survive after 
CCSN. Therefore, the birthrate and number in 
case 2 are about 1/3 of these in case 1, and 
there is not PLMXBs with RG (Wind) donors. 
In cases 1 and 3, different algorithms of CE are 
used. Usually, binary orbit after CE shortens up to 
~1% of initial one under a- algorithm assumption, 
while it approximately remains unchanged under 
7-algorithm assumption. Many binaries avoid to 
merge when they are undergoing CE evolution in 
case 3. There are more PLMXBs in case 3 than 
those in case 1. 
20 to 10 km s-i 
cr^ = 20 km s 

than those with ajj = 10 km s^^ for the progen 



Parameter is decreased from 

in case 5. AIC and EIC with 

^ produce wider orbital periods 
1 



itors of PLMXBs although some binary systems 
can be disturbed. Wider orbital periods provide 
enough separations so that the secondaries can 
evolve to RG, and can survive after CE evolution. 
Therefore, PLMXBs via AIC and EIC in case 5 are 
lower than those in case 1. We also carried out a 
test in which crj — 190 km s~^, and found that 
the number and birthrate of PLMXBs via AIC 
and EIC greatly decrease. Compared cases 6 and 
7 with case 1, parameter has a weak effect on 
PLMXBs populations. Therefore, from this later, 
we only discuss cases 1, 2, 3 and 5 to illustrate the 

Table 1: Parameters of the models for PLMXBs' populatiS fe^s of input parameters on PLMXBs, respec- 

tively. 

Eq.([T|) The number of PLMXBs mainly depends on the 
Eq.([T|) masses and the evolution phase of donors. Figure[T] 
Eq.((T]) shows the distribution of donors' masses. Majority 
Eq.([T]) of MS donors have lower masses than 1 Mq, while 
Eq. Ill) most of RG (Wind) donors have masses than 1 Mq 
= 1.2and the durations of RGs are very shorter than 
qc = 2.0those of MSs. Therefore, majority of PLMXBs are 
NS + MS systems, and PLMXBs with RG (Wind) 



cases 
case 1 
case 2 
case 3 
case 4 
case 5 
case 6 
case 7 



(km s ^) 


CE 


aj:(km s-^) 


190 


ttcoAcc = 1.0 


20 


400 


accAcc = 1.0 


20 


190 


7 = 1.5 


20 


400 


7 = 1.5 


20 


190 


ttceAcc = 1.0 


10 


190 


acoAce = 1.0 


20 


190 


accAce = 1.0 


20 
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Table 2: Different models of PLMXBs' population. The first column gives the model number according to Table 
[TJ Columns 2 to 11 give the birthrates and the numbers of PMLXBs with different kinds of donors, respectively. 
Total birthrate and number are showed in columns 12 and 13, respectively. NS + MS means that accreting 
NS has a MS donor, NS + RG (Roche) and NS + RG (Wind) mean that accreting NS is fed via Roche lobe 
overflow and stellar wind from a RG donor, NS + WD represents that accreting NS has a WD donor, and NS 
+ He represents that accreting NS has a naked He donor, respectively. 

"Cases NS+MS NS+RG NS+RG NS+He NS+WD Total 

Roche Wind 



BIR 




NUM 


BIR 




NUM 


BIR 




NUM 


BIR 




NUM 


BIR 




NUM 


BIR 




NUM 


(1) (2) 




(3) 


(4) 




(5) 


(6) 




(7) 


(8) 




(9) 


(10) 




(11) 


(12) 




(13) 


CCSN 




































case 1 4.1x10" 


-5 


68000 


1.1 X 10" 


-5 


830 


4.6 X 10" 


-6 


250 


3.7 X 10 


-5 


110 


1.1 X 10" 


-5 


50 


1.1 X 10" 


-4 


70000 


case 2 1.5x10" 


-5 


23000 7.4 X 10- 


-6 


460 


2.7 X 10" 


-7 


< 10 4.2 X 10" 


-6 


20 


3.7 X 10" 


-6 


30 


2.5 X 10" 


-5 


24000 


case 3 7.0x10" 


-5 


1000004.1 X 10" 


-5 


2000 6.4 X 10" 


-6 


140 


1.7 X 10" 


-5 


50 


9.3 X 10" 


-6 


50 


1.3 X 10" 


-4 


110000 


case 4 5.2x10" 


~5 


78000 3.0 X 10- 


-5 


1200 2.0 X 10" 


-6 


30 


3.6 X 10 


-6 


10 


1.6 X 10" 


'6 


10 


8.2 X 10" 


-5 


80000 


case 5 4.1x10" 


~5 


68000 


1.1 X 10" 


-5 


830 


4.6 X 10" 


-6 


250 


3.7 X 10 


-5 


110 


1.1 X 10" 


-5 


50 


1.1 X 10" 


-4 


70000 


case 6 4.1x10" 


~5 


68000 


1.2 X 10- 


-5 


840 


4.6 X 10" 


-6 


250 


3.6 X 10" 


-5 


90 


9.5 X 10- 


-6 


40 


1.4 X 10" 


-4 


70000 


case 7 4.0x10" 


-5 


66000 


1.5 X 10- 


-5 


850 


4.6 X 10" 


-6 


250 


3.4 X 10" 


-5 


70 


6.1 X 10" 


-6 


30 


1.1 X 10" 


-4 


68000 


AIC 




































case 1 2.1x10" 


-6 


180 


1.1 X 10" 


-5 


300 








3.4 X 10" 


-4 


1700 2.0 X 10" 


-4 


470 


4.7 X 10" 


-4 


2700 


case 2 2.1x10" 


-6 


180 


1.1 X 10" 


-5 


300 








3.4 X 10" 


-4 


1700 1.9 X 10- 


-4 


450 


4.6 X 10" 


-4 


2600 


case 3 4.2x10" 


-6 


470 


2.0 X 10" 


-5 


340 








6.2 X 10" 


-5 


280 


6.9 X 10" 


-5 


280 


1.3 X 10" 


-4 


1400 


case 4 4.2x10" 


-6 


460 


2.0 X 10" 


-5 


340 








6.2 X 10" 


-5 


280 


6.9 X 10" 


-5 


300 


1.3 X 10" 


-4 


1400 


case 5 2.1x10" 


-6 


220 


1.1 X 10" 


-5 


280 








3.4 X 10" 


-4 


1700 2.0 X 10" 


-4 


490 


4.7 X 10" 


-4 


2700 


case 6 2.7x10" 


-6 


180 


2.5 X 10" 


-5 


320 








2.8 X 10 


-4 


1300 1.5 X 10" 


-4 


420 


3.9 X 10" 


-4 


2200 


case 7 2.7x10" 


-6 


180 


9.4 X 10" 


'5 


500 








1.3 X 10" 


-4 


530 


8.6 X 10" 


^4 


300 


2.8 X 10" 


-4 


1500 


EIC 




































case 1 1.5x10" 


-6 


1100 


2.9 X 10" 


-5 


140 


4.2 X 10" 


-6 


< 10 1.1 X 10 


-4 


440 


2.8 X 10" 


-5 


220 


1.4 X 10" 


-4 


1900 


case 2 1.5x10" 


-6 


1100 


2.9 X 10" 


-5 


140 


4.2 X 10" 


-6 


< 10 1.1 X 10" 


-4 


440 


2.8 X 10" 


-5 


220 


1.4 X 10" 


-4 


1900 


case 3 3.5x10" 


-7 


360 


2.2 X 10" 


-5 


40 


7.1 X 10" 


-6 


< 10 6.6 X 10" 


-5 


460 


4.9 X 10" 


-5 


350 


1.2 X 10" 


-4 


1200 


case 4 3.5x10" 


-7 


360 


2.2 X 10" 


-5 


40 


7.1 X 10" 


-6 


< 10 6.6 X 10 


-5 


460 


4.9 X 10" 


-5 


350 


1.2 X 10" 


-4 


1200 


case 5 7.5x10" 


-7 


850 


2.6 X 10" 


-5 


90 


4.5 X 10" 


-6 


< 10 9.4 X 10" 


-5 


450 


3.0 X 10" 


-5 


210 


1.4 X 10" 


-4 


1600 


case 6 1.5x10" 


-6 


1100 


2.9 X 10" 


-5 


140 


4.2 X 10" 


-6 


< 10 1.0 X 10" 


-4 


440 


2.7 X 10" 


-5 


210 


1.7 X 10" 


-4 


1900 


case 7 1.5x10" 


-6 


1100 


2.9 X 10" 


-5 


150 


4.1 X 10" 


-6 


< 10 


8.9 X 10" 


-5 


310 


2.3 X 10" 


-5 


190 


1.8 X 10" 


-4 


1800 
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donors are very rare. 

In PLMXBs with WD donors, matter transfer is 

driven by K ravitational radiatio n. Delove fc BildstenI 



(|2003l ) and iBildsten fc Delovel (I200J) showed the 
relation among WD masses, orbital periods and 
mass-transfer rates in ultra-compact X-ray bina- 
ries with WD donors. According to their results, 
the higher WD donors' masses are, the higher 
mass-transfer rates are. The duration of PLMXBs 
with massive WDs is very short. As Figure [T] 
shows, there are two peaks in the distributions of 
WDs' masses. The left peak is at ~ 0.03 Mq, 
and the right peak is at ^ 0.09 Mq. The former 
mainly comes from the PLMXBs which will trans- 
late from persistent to transient state because the 
donors in these PLMXBs have low mass-loss rates 
which results in long durations. The later mainly 
results from the PLMXBs who have undergone the 
AIC. Compared with the PLMXBs around the left 
peak, these PLMXBs around the right peak have 
short orbital periods and high X-ray luminosi- 
ties. There are ~ 600 — 900 PLMXBs with WD 
donors. However, the duration of LMXBs with 
WD donors whose masses are lower than those in 
PLMXBs are very long (~ 10^ yr) because of a low 
mass-transfer rate (~ 1O~^^M0 yr""'^). Therefore, 
most of LMXBs with WD donors are transient. 

Not like WDs, naked He stars are convective. 
The mass transfer in PLMXBs with He donors 
is driven by gravitational radiation and mag- 
netic braking. In general, the later is dominate d 
(|Hurlev et all l2002t IPostnov fc Kuranovl l2005f ). 
and drives a mass-transfer rate of ~ 10"'' Mq yr"-'^ 
in our work. Majority of He donors' masses are 
between ~ 0.3 — 2.0 Mq (Figure [J). There- 
fore, there are several thousand PLMXBs with 
He donors in the Galaxy. According to Table [H 
the number of PLMXBs with He donors is much 
larger than PLMXBs with WD d onors. From the 



prope rties of type I X-ray bursts, lin't Zand et al 



(j2005l ) suggested that in most ultra-compact X- 
ray binaries the matter accumulated on NSs is 
helium. This is consistent with our results al- 
though we do not discuss transient LMXBs. 
In the PLMXBs plotted in Figure [21 there are 
two ultra-compact X-ray binaries (4U 1626-67 
and 4 U 0614-1-09) which ha ve very evolved He 
donors! Nelemans et al. 2010t) . 



3.2. Properties of PLMXBs with Different 
Donors 

As SJT] mentions, the donors of LMXBs basically 
determine the orbital periods and mass-transfer 
rates which give the X-ray luminosity. Figure [2] 
shows the distributions of the orbital periods and 
the X-ray luminosities or the mass-accretion rate 
Mns- PLMXBs with different donors have differ- 
ent positions in Figure [H 

Revnivtsev et al. ( 201 if ) investigated the bright- 



est Galactic PLMXBs which are plotted in Figure 
El They concluded that the majority of PLMXBs 
with X-ray luminosities below 2 x 10'^'' erg s~^ 
have unevolved MS, while PLMXBs with higher 
X-ray luminosity predominantly harbor giant 
donors. In the panels of NS + MS in Figure 121 ~ 
90% of PLMXBs with MS donors lie in the shallow 
region with a X-ray luminosity of ^ lO'^^ erg s~^. 
The irradiation of LMXBs can drive mass tra nsfer 
(|Podsiadlowskilll99ll: [Biining fc Ritteil |2004| ). In 



this work, we do not consider the effect of irradi- 
ation. Therefore, compared with the luminosities 
of known PLMXBs, we may underestimate the 
luminosity of PLMXBs with MS donors. In the 
panels of NS + RG (Roche) in Figure |21 our sam- 
ple covers the positions of Cyg X-2 and GX 13+1 
which have long orbital periods. These PLMXBs 
can only be explained by N S -|- RG (Roche) 



Orosz fc Kuulker j 



system s in our simulations. 

()l999l) gave good measurements for Cyg X-2 
and have derived the donor's mass around 0.6 

Mq. 



Podsiadlowski fc Rappaport ( 2000l ) sug 



gested that the donor in Cyg X-2 has a mass 
of around 0.5 Mq with a non-degenerated helium 
core and is burning hydrogen in a shell. That is, 
the donor in Cyg X-2 is a sub-giant which has un- 
dergone t he violent mass loss. This is consistent 
with ours. iBandvopadhvav et al.l (|l999f ) gave that 
the spectrum of donor in GX 13-1-1 clearly shows 
the features of the K5 HI giant. Figure [31 gives the 
distributions of donors' luminosities and mass- 
transfer rates determine the X-ray luminosity. As 
the panels of NS -I- RG (Roche) and NS -I- RG 
(Wind) in Figure |3l show, the donors' luminosities 
in PLMXBs with RG (Roche) donors are much 
lower than those in PLMXBs with RG (Wind), 
while the X-ray luminosities in the former are 
much higher than those in the later. Therefore, it 
is more difficult to observe the donors' luminosi- 
ties in PLMXBs with RG (Roche) donors. We 
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suggested that donors in Cyg X-2 and GX 13+1 
are giants which fill up Roche lobes. The orbital 
periods and X-ray luminosities of 2 known sym- 
biotic X-ray binaries (GX 1+4 and 4U 1700+24) 
are measured, which is plotted by triangles in 
Figure [2] NS + RG (Wind) systems can cover 
4U 1700+24 very well, but our models can not 
explain GX 1+4. A detailed investigation of sym- 
biotic X-ray binaries is being carried out (Lii et 
al. in preparation). 

In our simulations, the orbital periods of 
PLMXBs with WD or He donors are shorter than 
80 minutes, and they are ultra-compact X-ray 
binaries. About 10% — 34% of WDs in NS 
+ WD systems are He WDs, and 66% — 90% 
are CO WDs. NS + ONeMg WD systems are 
negligible. As Figure [2] shows, PLMXBs with 
WD donors are agree with observations, while 
PLMXBs with He donors have X-ray luminosi- 
ties higher than those of observ ed ult ra-compact 
X-ray binaries. Nelemans et al\ ( 2010t ) suggested 
that two ultra-compact X-ray binaries 4U 1626-67 
and 4U 0614+09 have very evolved He donors, 
and their X-ray luminosities are 3.2 x 10'^^ erg s~^ 
and 3.4 x 10^^ erg s^^. In our model, the mass 
transfer in PLMXBs with He donors is driven by 
the magnetic braking, which produces the X-ray 
luminosity of ^ 10'^* erg s~^. We may overesti- 
mate the work efficiency of the magnetic braking 
driving mass transfer. As panels of NS + WD 
and NS + He in Figure [3] show, the donors' lu- 
minosities in PLMXBs with He donors are much 
higher than those in PLMXBs with WD donors. 
This difference may be a way via which we can 
distinguish He donors from WD donors. 

Figure S] gives the distribution of the X-ray 
luminosities (mass-accretion rates) of NSs in 
PLMXBs with different donors. The X-ray lu- 
minosities between ~ 6 x 10^^ — 6 x 10'^'^ erg 
s"""^ mainly come from PLMXBs with RG (Wind). 
The X-ray luminosities between ^ 10^^ — 10^^ 
erg s~i originate from PLMXBs with MS, RG 
(Roche) and WD donors, in which PLMXBs 
with MS are dominated. The X-ray luminosi- 
ties between ~ 10^^ — 10^^ erg s~^ mainly orig- 
inate from PLMXBs with He donors in which 
the mass transfer is driv e n by magnetic brak- 



the accretion of matter onto a NS with magnetic 
stellar wind, which agrees with our results. How- 
ever, the observed X-ray luminosity function ( up 
to - 2 x 10^^ erg s^^ ) of LMXBs can be explained 
by PLMXBs with MS and RG (Roche) donors. In 
BSE code, gravitational radiation is only efficient 
for binaries with orbital periods less than 3 hours. 
In order to explain a gap of CVs between 2 and 3 
hr in the otherwise smooth period-mass distribu- 
tion, BSE code does not apply magnetic braking 
when the primary is a fully convective MS star 
whose mass is lower than 0.35Mq. There is not 
the magnetic braking in NS + MS systems if MS's 
mass is larger than 1.2M0 because it has no con- 
vective envelope. In our work, less than 30% of 
PLMXBs with MS donors have a donor whose 
mass is between 0.35Mo and 1.2M0. Therefore, 
the mass transfer in PLMXBs with X-ray luminos- 
ity function ( up to ^ 2 x 10'^^ erg s~^ ) is mainly 
driven by stellar evolution or magnetic braking, 
but not gravitational radiation. 

In the bottom panels of Figure HI we give 
the distributions of total X-ray luminosities which 
are different from the power-law X-ray luminos- 
ity funct i on (sh owed by dashed line) observed by 
iGilfanovl (|2004l) . There are two main reasons: 
i)Our work does not include the transient X-ray 
luminosity. The transient systems in outburst 
reach high (close to Eddington) X-ray luminosi- 



ing. Postnov fc Kuranovl ( 2005 ) suggested that 



observed X-ray luminosity function (> 2 x 10^^ 
erg s~^ ) of LMXBs can generally be explained by 



ties. iBelczvnski et a l. (2008) suggested that the 
X-ray luminosities in outburst are between ~ 10^^ 
and 10'^^ erg s~^. ii)We overestimate the X-ray lu- 
minosities of PLMXBs with He donors. Of course, 
we must obtain total population of LMXBs in- 
cluding transient and persistent systems with NSs 
and black hole before we compare X-ray luminos- 
ity function with observational one. In the future, 
we will try doing it. 

4. Conclusion 

We perform a detailed study of donors in 
PLMXBs, employing the population synthesis ap- 
proach to the evolution of binaries. We estimate 
that there are - 29000 ~ 110000 PLMXBs in 
the Galaxy, and their birthrates are ~ 3.4 — 7.2 
xl0~^ yr~^. PLMXBs with different donors have 
different formation channels. Our numerical sim- 
ulation shows that more than 90% of PLMXBs 
have MS donors, and PLMXBs with RG (Wind) 
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Fig. 2. Distributions of the orbital periods vs. the X-ray luminosities or the mass-accretion rate Mns 

in PLMXBs. The gradations of gray-scale correspond to the regions where the number density of systems 



is, respectively, within 1 - 1/2, 1/2 - 1/4, 1/4 - 1/8, 1/8 - of the maximum of 



-. Number in 



every panel is normalized to 1 . Circles represents th e brightest PL MXBS from iRevnivtsev et al.l (|201lh , and 
triangles represent the symbiotic X-ray binaries in Masetti et al. (j2007l) . Solids represent the relation for 
hydrogen-rich disk of Eq.([5l), and dashed, dotted and dash-dotted lines represent the relations for He-rich, 
C-rich and 0-rich disk of Eq.([6]), respectively. 
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donors arc negligible. Most of NSs in PLMXBs 
with hydrogen- rich donors have undergone CCSN, 
while more than 90% of NSs in PLMXBs with He 
donors or WD donors have undergone AIC and 
EIC. 

In PLMXBs with MS donors, the orbital pe- 
riods are between ~ 1 hour and 100 hours, and 
the mass transfer is driven by donor evolution 
or magnetic braking. Our population synthesis 
code shows that their X-ray luminosities mainly 
are around ~ 10^^ erg s~^. Similarly, in PLMXBs 
with RG donors via Roche lobe overflow(Roche), 
the mass transfer is driven by donor evolution, 
but orbital periods are between ~ 10 hours and 
1000 hours. Their X-ray luminosities are around 
- 10'^^ erg s-i. The 2 known LMXBs ( Cyg X-2 
and GX 13+1) can belong to PLMXBs with RG 
(Roche). PLMXBs with RG (Wind) donors have 
the longest orbital periods and low X-ray luminosi- 
ties. Their contributions to X-ray luminosities can 
be negligible. In PLMXBs with He donors, the 
orbital periods are shorter than 80 minutes, and 
the mass transfer is mainly driven by magnetic 
braking. Results of our numerical simulations pre- 
dict that PLMXBs with X-ray luminosities around 
~ 10^^ erg s~^ mainly come from binaries with 
He donors, but their X-ray luminosities may be 
overestimated in our work. In PLMXBs with WD 
donors, the orbital periods are shorter than 1 hour, 
and the mass transfer is mainly driven by gravita- 
tional radiation. According to results of our popu- 
lation synthesis code, their X-ray luminosities are 
between ~ 6 x 10^^ — 10^^ erg s~^, and most of 
LMXBs with WD donors are transient. 

In this work we do not consider transient 
LMXBs with different donors. In hirtlier work, we 
will investigate donors in persistent and transient 
LMXBs, and discuss X-ray luminosity function of 
LMXBs. 
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Donors of Persistent Neutron-star Low-mass X-ray Binaries 



Chunhua Zhu^'^, Guoliang Lii^, Zhaojun Wang^, Na Wang^ 
^ chunhuazhuO s ina . cn ; na . wangSxao .ac.cn 

ABSTRACT 

Properties of X-ray luminosities in low-mass X-ray binaries (LMXBs) mainly depend 
on donors. We have carried out a detailed study of donors in persistent neutron-star 
LMXBs (PLMXBs) by means of a population synthesis code. PLMXBs with different 
donors have different formation channels. Our numerical simulations show that more 
than 90% of PLMXBs have main sequence (MS) donors, and PLMXBs with red giant 
(RG) donors via stellar wind (Wind) are negligible. In our model, most of neutron stars 
(NSs) in PLMXBs with hydrogen-rich donors form via core-collapse supernovae, while 
more than 90% of NSs in PLMXBs with naked helium star (He) donors or white dwarf 
(WD) donors form via an evolution-induced collapse via helium star (1.4 < Mus/Mq < 
2.5) or an accretion-induced collapses for an accreting ONeMg WD. 

PLMXBs with different donors have different properties. In PLMXBs with MS 
donors, the orbital periods are between ~ 1 hour and 100 hours, and the mass transfer 
is driven by donor evolution or magnetic braking. Our population synthesis code shows 
that their X-ray luminosities mainly are around ~ 10^^ erg s^^. Similarly, in PLMXBs 
with RG donors via Roche lobe overflow(Roche), the mass transfer is driven by donor 
evolution, but orbital periods are between ~ 10 hours and 1000 hours. Their X-ray 
luminosities are around ~ lO^^ erg s'^. The 2 known LMXBs ( Cyg X-2 and GX 
13+1) can belong to PLMXBs with RG (Roche). PLMXBs with RG (Wind) donors 
have the longest orbital periods and low X-ray luminosities (~ 10^^ erg s~^). Their 
contributions to X-ray luminosities can be negligible. In PLMXBs with He donors, the 
orbital periods are shorter than 80 minutes, and the mass transfer is mainly driven 
by magnetic braking. Results of our numerical simulations predict that PLMXBs with 
X-ray luminosities around ~ 10'^® erg mainly come from binaries with He donors. 
In PLMXBs with WD donors, the orbital periods are shorter than 1 hour, and the 
mass transfer is mainly driven by gravitational radiation. According to results of our 
population synthesis code, their X-ray luminosities are between ~ 6 x 10^^ — 10^^ erg 
s~^, and most of LMXBs with WD donors are transient. 

Subject headings: h 
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inaries: close — pulsar: general — stars: neutron — X-ray: stars 

1. Introduction 

Low-mass X-ray binaries (LMXBs) were discovered nearly 50 yr ago, and there are now ^ 

200 known in the Galaxy (?). LMXB is a mass-transferring binary system with a compact object 
accretor(a black hole or a neutron star (NS)), and a low-mass (< IMq) donor. The X-ray luminosity 
function is an important characteristic of the LMXBs, and has become a key tool for studying 
LMXBs. 

Using results of Chandra observations of old stellar systems in 11 nearby galaxies of various 
morphological types and the census of LMXBs in the Galaxy, ? suggested that the total number 
of LMXBs and their combined luminosity are proportional to the stellar mass of the host galaxy. 
? suggested that the flattening of the LMXBs luminosity function at lower than 2 x 10^^ erg s~-^ 
might correspond to the transition from the magnetic stellar wind braking to the gravitational 
wave braking mechanism. ? suggested that LMXBs with X-ray luminosities below 2 x 10'^^ have 
unevolved secondary companions while systems with higher X-ray luminosity predominantly harbor 
giant donors. 

In theoretical work, people usually assume that the X-ray luminosity is directly proportional 
to the mass-accretion rate of compact stars in LMXBs. The mass-accretion rate depends on the 
orbital period, donor's evolution, and angular momentum loss. Donor is usually a main sequence 
(MS) or a white dwarf (WD), and it may also be a red giant (RG) or a naked helium star (He). The 
properties of LMXBs are closely related to their donors. If donors are poor-hydrogen stars (WDs 
or Hes), LMXBs usually are ultra-compact X-ray binaries whose orbital periods are shorter than 
80 minutes (?). If donors are RGs, LMXBs are called as symbiotic X-ray binaries(?). Up to now, 
there are about a dozen of ultra-compact X-ray binaries and about 10 symbiotic X-ray binaries. 

Different types of stars have different properties and evolutions. Therefore, the donors in 
LMXBs determine orbital period and angular momentum loss, which directly affects the mass- 
transfer rate from the donors to NSs. Obviously, the donors really determine the luminosity of 
LMXBs. Simultaneously, the observational luminosity function of LMXBs provides important 
constrains for our simulating donors' evolution. 

On observations, LMXBs are divided into transient and persistent sources. It is difficult to 

estimate the X-ray luminosities of transient LMXBs during quiescent state and outburst state. 
In this paper, we focus on persistent LMXBs (PLMXBs) with accreting NS and different donors, 
and investigate their properties and contributions to total X-ray luminosity. In § 2 we present 
our assumptions and describe some details of the modeling algorithm. In § 3 we discuss the main 
results. In § 4 the main conclusions are given. 
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2. Models 

For the simulation of binary evolution, we use rapid binary star evolution code BSE (?) with 
updates by ?. BSE code calculates the orbital changes of binary systems via mass variations, 
gravitational radiation and magnetic braking. Details are in §2 of ?. 

Forming channels of LMXBs are important questions in X-ray astronomy. There are many 
literatures to investigate them (e.g., ????). The difficulty in a theoretical investigation involves two 
problems: (i)keeping the binary bound when the massive progenitor of NS explodes in a supernova 
event; (ii)a common envelope (CE) phase. They have great effects on LMXBs. The following two 
subsections give descriptions. 



2.1. Common Envelope Evolution 

In binary system, due to orbital angular momentum loss or stellar expansion, a star can 
overflow its Roche lobe. If the mass ratio of the components {q = MdonorZ-^accretor) at the onset 
of Roche lobe overflow is larger than a certain critical value Qc, the mass transfer is dynamically 
unstable and results in the formation of a CE. The issue of the criterion for dynamically unstable 
Roche lobe overflow, qc, is still open. Based on the polytropic models, ? gave qc for red giants by 

^^ = °-^^^ + 3x(l-i/M,_)' 
where Mc is the core mass of donor. However, this qc is obtained under conservative Roche lobe 
overflow. ?? showed that qc depends heavily on the assumed mass-transfer efficiency. They found 
that qc almost linearly increases with the amount of the mass and momentum lost during mass 
transfer. In ? the critical mass ratio qc is between 1.1 and 1.3. In this work, we take qc =Eq.([T]), 
qc = 1.2 and qc = 2.0 in different simulations, respectively. 

Although many efforts have been devoted to understanding the evolution of CE (e.g., ???), 
the knowledge about it is still poor. It is generally assumed that the orbital energy of the binary 
is used to expel the envelope of the donor with an efficiency ace, which is called as a-algorithm. 
? suggested to describe the CE evolution by an algorithm based on the equation for the system 
orbital angular momentum balance which implicitly assumes the conservation of energy(?), which 
is called as 7-algorithm. Following ?, for CE evolution in different simulations we use ace^ce = 1-0 
in a-algorithm and 7 = 1.5 in 7-algorithm, respectively. Here Ace is a structure parameter which 
depends on the evolutionary stage of the donor. 



2.2. Formation Channels of Neutron Stars and Kick Velocity 



In X-ray binaries, NSs can be formed via three channels(e.g., ??): (i) Core-collapse supernovae 
(CCSN) for a star with main-sequence mass M/Mq > 11; (ii)Evolution induced collapse (EIC) of 
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a helium star with a mass between 1.4 and 2.5Mq in which the cohapse is triggered by electron 
capture on ^'^Ne and -^"^Mg (?); (iii) Accretion-induced collapses (AIC) for an accreting ONeMg 
WD whose mass reaches the Chandrasekhar limit. Response of accreting ONeMg WD is treated 
in the same way as the evolution of CO WD (see details in ?). 

Nascent NS receives additional velocity ( "kick" ) due to some still unclear process that disrupts 
spherical symmetry during the collapse or later Dichotomous nature of kicks which was suggested 
quite early by ?. Observationally, the kick is not well constrained due to numerous selection effects. 
Currently, high kicks (~ 100km s^^) are associated with NS originating from CCSN, while low 
kicks (~ 10km s"^) with NS born in EIC and AIC (?). 

We apply to core-collapse NS Maxwellian distribution of kick velocity t^k 




o-k = 190 and 400 km s"^ for CCSN, while = 20 and 10 km s"^ for EIC and AIC in different 
simulations. 



2.3. X-ray Luminosity 

The X-ray luminosity of the accreting NS can be approximated by 

Lboi = r/MNSC^ ^ 5.7 X 1035erg s~H^){j^^=J^^) (3) 

where rj ~ 0.1 is the efficiency of accretion onto the NS and Mns is the mass-accretion rate of 
the NS. Super-Eddington accretion rates may be important in the formation of low- mass X-ray 
binaries and millisecond pulsars(?). We assume that Mns = min(MNS) 'yEdd ^ -^Edd)) where MEdd 
is Eddington limit given by 

MEdd = 2.08 X 10-3(1 + X)-'R^sMQyi-\ (4) 

Here, X is the hydrogen mass fraction, rysdd is the factor to allow super-Eddington luminosities, 
taken to be 5 (??). To transform the bolometric luminosity into the X-ray luminosity, a bolometric 
correction factor ryboi is introduced by Lx = f?bol-^bol- Following ?, we take r/boi = 0.55. 

Roche overflow-fed systems are subject to a thermal disk instability and may appear either 
as persistent or transient X-ray sources depending on the mass transfer rate. A system becomes a 
transient X-ray source when the mass-transfer rate falls below a certain critical value. Merit. For 
hydrogen-rich disks (The donors are MSs or RGs ), we use the work of ?. Applying to Eq.Q, Merit 
for hydrogen-rich disk is given by 

Merit = 1.8 X 10i5pi,{]^ g/s, (5) 
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where Porb is orbital period in hours. For disks with heavier elements, we use the work of ?: 

( 5.9 X lOi^MNs'^^i^i^^ag-f g/s, He rich 

Merit ={ 1.2 X 10i6M-0-7X''«o:f g/s, Crich (6) 
[ 5.0 X 1016m-°-68 ^2.05^0.45 g/g^ o rich 

where Rd is a maximum disk radius (2/3 of accretor Roche lobe radius) in 10^'^ cm, ao.i = a/0.1 
in which a = 0.1 is a viscosity parameter. 

If Mns > -M^crit or wind-fed accretion, the system is a PLMXB whose X-ray luminosity is 
determined by Eq. ([3]). If Mns < -M'crit in Roche overflow accretion, the system is a transient 
source. In this work we focus on PLMXBs. 

3. Results 

We use Monte Carlo method to simulate the initial binaries. For initial mass function, mass- 
ratios, and separations of components in binary systems, we adopt the distributions used by us in 
??. We assume that all binaries have initially circular orbits. After a supernova, new parameters 
of the orbit are derived using standard formulae, (e. g., ?). Table [T] lists all cases considered in 
the present work. Our model is normalized to formation of one binary with Mi > 0.8Mq per year 
(?). We use 1 X 10^ binary systems in our Monte-Carlo simulations. 

In this work, a binary is considered as PLMXB if it satisfies the following conditions: 

(i) Binary includes an NS and its companion's mass is lower than 6 Mq; 

(ii) Binary orbital period and the mass-accretion rate of NS satisfy Eqs.([5]) and ([6]). 
Here, we call both low- and intermediate-mass X-ray binaries as LMXBs. 

Table 1: Parameters of the models for PLMXBs' populations. 



cases 


crk(km s ^) 


CE 


alikm s~^) 




case 1 


190 


"ceAcc = 1.0 


20 


Eq.© 


case 2 


400 


"ceAcc = 1-0 


20 


Eq.© 


case 3 


190 


7 = 1.5 


20 


Eq.© 


case 4 


400 


7 = 1.5 


20 


Eq.© 


case 5 


190 


"ceAcc = 1-0 


10 


Eq.© 


case 6 


190 


"ceAce = 1-0 


20 


Qc = 1.2 


case 7 


190 


OccAcc = 1-0 


20 


Qc = 2.0 
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Table 2: Different models of PLMXBs' population. The first column gives the model number according 
to Table [TJ Columns 2 to 11 give the birthrates and the numbers of PMLXBs with different kinds of 
donors, respectively. Total birthrate and number are showed in columns 12 and 13, respectively. NS 
+ MS means that accreting NS has a MS donor, NS + RG (Roche) and NS + RG (Wind) mean that 
accreting NS is fed via Roche lobe overflow and stellar wind from a RG donor, NS + WD represents 
that accreting NS has a WD donor, and NS + He represents that accreting NS has a naked He donor, 
respectively. 

Cases NS+MS NS+RG NS+RG NS+He NS+WD Total 

Roche Wind 

BIR NUM BIR NUM BIR NUM BIR NUM BIR NUM BIR NUM 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 
CCSN 



case 1 4.1x10 


_5 


68000 


1.1 X 10" 


_5 


830 


4.6 X 10" 


_g 


250 


3.7 X 10" 


_5 


110 


1.1 X 10" 


_5 


50 


1.1 X 10" 


-4 


70000 


case 2 1.5 x lU 


-5 


23000 7.4 X 10" 


-6 


460 


2.7 X 10" 


-7 


< 10 4.2 X 10" 


-6 


20 


3.7 X 10" 


-6 


30 


2.5 X 10 


-5 


24000 


case 3 7.0x10' 


-5 


100000 4.1 X 10" 


-5 


2000 6.4 X 10" 


-6 


140 


1.7 X 10" 


-5 


50 


9.3 X 10" 


-6 


50 


1.3 X 10" 


-4 


110000 


case 4 5.2x10" 


-5 


78000 3.0 X 10" 


-5 


1200 2.0 X 10" 


-6 


30 


3.6 X 10" 


-6 


10 


1.6 X 10" 


-6 


10 


8.2 X 10" 


-5 


80000 


case 5 4.1x10" 


-5 


68000 


1.1 X 10" 


-5 


830 


4.6 


X 10" 


-6 


250 


3.7 X 10" 


-5 


110 


1.1 X 10" 


-5 


50 


1.1 X 10" 


-4 


70000 


case 6 4.1x10" 


-5 


68000 


1.2 X 10" 


-5 


840 


4.6 X 10" 


-6 


250 


3.6 X 10" 


-5 


90 


9.5 X 10" 


-6 


40 


1.4 X 10" 


-4 


70000 


case 7 4.0x10" 


-5 


66000 


1.5 X 10" 


-5 


850 


4.6 X 10" 


-6 


250 


3.4 X 10" 


-5 


70 


6.1 X 10" 


-6 


30 


1.1 X 10" 


-4 


68000 


AIC 






































case 1 2.1x10" 


-6 


180 


1.1 X 10" 


-5 


300 










3.4 X 10" 


-4 


1700 2.0 X 10" 


-4 


470 


4.7 X 10" 


-4 


2700 


case 2 2.1x10" 


-6 


180 


1.1 X 10" 


-5 


300 










3.4 X 10" 


-4 


1700 1.9 X 10" 


-4 


450 


4.6 X 10" 


-4 


2600 


case 3 4.2x10" 


-6 


470 


2.0 X 10" 


-5 


340 










6.2 X 10" 


-5 


280 


6.9 X 10" 


-5 


280 


1.3 X 10" 


-4 


1400 


case 4 4.2x10" 


-6 


460 


2.0 X 10" 


-5 


340 










6.2 X 10" 


-5 


280 


6.9 X 10" 


-5 


300 


1.3 X 10" 


-4 


1400 


case 5 2.1x10" 


-6 


220 


1.1 X 10" 


-5 


280 










3.4 X 10" 


-4 


1700 2.0 X 10" 


-4 


490 


4.7 X 10" 


~4 


2700 


case 6 2.7x10" 


-6 


180 


2.5 X 10" 


-5 


320 










2.8 X 10" 


-4 


1300 1.5 X 10" 


-4 


420 


3.9 X 10" 


-4 


2200 


case 7 2.7x10" 


-6 


180 


9.4 X 10" 


-5 


500 










1.3 X 10" 


-4 


530 


8.6 X 10" 


-4 


300 


2.8 X 10" 


-4 


1500 


EIC 






































case 1 1.5x10" 


-6 


1100 


2.9 X 10" 


-5 


140 


4.2 


X 10" 


-6 


< 10 1.1 X 10" 


-4 


440 


2.8 X 10" 


-5 


220 


1.4 X 10" 


-4 


1900 


case 2 1.5x10" 


-6 


1100 


2.9 X 10" 


-5 


140 


4.2 


X 10" 


-6 


< 10 1.1 X 10" 


-4 


440 


2.8 X 10" 


-5 


220 


1.4 X 10" 


-4 


1900 


case 3 3.5 x 10" 


-7 


360 


2.2 X 10" 


-5 


40 


7.1 


X 10" 


-6 


< 10 6.6 X 10" 


-5 


460 


4.9 X 10" 


-5 


350 


1.2 X 10" 


-4 


1200 


case 4 3.5x10" 


-7 


360 


2.2 X 10" 


-5 


40 


7.1 


X 10" 


-6 


< 10 6.6 X 10" 


-5 


460 


4.9 X 10" 


-5 


350 


1.2 X 10" 


-4 


1200 


case 5 7.5x10" 


-7 


850 


2.6 X 10" 


-5 


90 


4.5 X 10" 


-6 


< 10 9.4 X 10" 


-5 


450 


3.0 X 10" 


~5 


210 


1.4 X 10" 


-4 


1600 


case 6 1.5x10" 


-6 


1100 


2.9 X 10" 


-5 


140 


4.2 X 10" 


-6 


< 10 1.0 X 10" 


-4 


440 


2.7 X 10" 


-5 


210 


1.7 X 10" 


-4 


1900 


case 7 1.5x10" 


-6 


1100 


2.9 X 10" 


-5 


150 


4.1 


X 10" 


-6 


< 10 8.9 X 10" 


-5 


310 


2.3 X 10" 


-5 


190 


1.8 X 10" 


-4 


1800 
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3.1. Birthrates and Numbers of PLMXBs' Populations 

In our simulations, there are ~ 29000 (case 2) — 110000 (case 3) PLMXBs in the Galaxy, and 
their birthrates are ~ 3.4 — 7.2 xlO~^ yr"^- However, the number of all observed LMXBs is less 
than 200 (?). ? investigated LMXBs via CCSN, and obtained birthrates for LMXBs of 10"^— 
10~^ yr~^ and 400 — 70000 LMXBs, a factor 10 — 1000 times higher than observed number. Many 
authors suspected that the mismatch between the observed number and theoretically predicted 
number could be related to irradiation effects (et al. ??). In this work, we do not consider LMXBs 
in low states driven by irradiation-driven limit cycles. However, we still encounter the known 
problem of overproduction of LMXBs. 

As Table [2] shows, NSs with different kinds of donors in PLMXBs have different formation chan- 
nels. More than 90% of PLMXBs have undergone CCSN, especially, for PLMXBs with hydrogen- 
rich donors. However, more than 90% of PLMXBs with He donors or WD donors have undergone 
AIC and EIC. In cases 1 and 2, parameter cTk is increased from 190 to 400 km s~^. The larger cjk 
is, the more difficultly binary survive after CCSN. Therefore, the birthrate and number in case 2 
are about 1/3 of these in case 1, and there is not PLMXBs with RG (Wind) donors. In cases 1 and 
3, different algorithms of CE are used. Usually, binary orbit after CE shortens up to ~1% of initial 
one under a-algorithm assumption, while it approximately remains unchanged under 7-algorithm 
assumption. Many binaries avoid to merge when they are undergoing CE evolution in case 3. There 
are more PLMXBs in case 3 than those in case 1. Parameter crjj is decreased from 20 to 10 km 
s~^ in case 5. AIC and EIC with = 20 km s~^ produce wider orbital periods than those with 
= 10 km for the progenitors of PLMXBs although some binary systems can be disturbed. 
Wider orbital periods provide enough separations so that the secondaries can evolve to RG, and 
can survive after CE evolution. Therefore, PLMXBs via AIC and EIC in case 5 are lower than 
those in case 1. We also carried out a test in which = 190 km s^^, and found that the number 
and birthrate of PLMXBs via AIC and EIC greatly decrease. Compared cases 6 and 7 with case 
1, parameter Qc has a weak effect on PLMXBs populations. Therefore, from this later, we only 
discuss cases 1, 2, 3 and 5 to illustrate the effects of input parameters on PLMXBs, respectively. 

The number of PLMXBs mainly depends on the masses and the evolution phase of donors. 
Figure [1] shows the distribution of donors' masses. Majority of MS donors have lower masses than 
1 Mq , while most of RG (Wind) donors have masses than 1 Mq and the durations of RGs are very 
shorter than those of MSs. Therefore, majority of PLMXBs are NS + MS systems, and PLMXBs 
with RG (Wind) donors are very rare. 

In PLMXBs with WD donors, matter transfer is driven by gravitational radiation. ? and ? 
showed the relation among WD masses, orbital periods and mass-transfer rates in ultra-compact 
X-ray binaries with WD donors. According to their results, the higher WD donors' masses are, 
the higher mass-transfer rates are. The duration of PLMXBs with massive WDs is very short. 
As Figure [U shows, there are two peaks in the distributions of WDs' masses. The left peak is at 
~ 0.03 Mq, and the right peak is at ~ 0.09 Mq. The former mainly comes from the PLMXBs 
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which will translate from persistent to transient state because the donors in these PLMXBs have 
low mass-loss rates which results in long durations. The later mainly results from the PLMXBs 
who have undergone the AIC. Compared with the PLMXBs around the left peak, these PLMXBs 
around the right peak have short orbital periods and high X-ray luminosities. There are ~ 600 — 
900 PLMXBs with WD donors. However, the duration of LMXBs with WD donors whose masses 
are lower than those in PLMXBs are very long (~ 10^ yr) because of a low mass-transfer rate 
(~ IO^^^Mq yr-i). Therefore, most of LMXBs with WD donors are transient. 

Not like WDs, naked He stars are convective. The mass transfer in PLMXBs with He donors 
is driven by gravitational radiation and magnetic braking. In general, the later is dominated (??), 
and drives a mass-transfer rate of ^ 10~'^Mq yr~^ in our work. Majority of He donors' masses 
are between ~ 0.3 — 2.0 Mq (Figured]). Therefore, there are several thousand PLMXBs with He 
donors in the Galaxy. According to Tabled the number of PLMXBs with He donors is much larger 
than PLMXBs with WD donors. Prom the properties of type I X-ray bursts, ? suggested that in 
most ultra-compact X-ray binaries the matter accumulated on NSs is helium. This is consistent 
with our results although we do not discuss transient LMXBs. In the PLMXBs plotted in Figure El 
there are two ultra-compact X-ray binaries (4U 1626-67 and 4U 0614+09) which have very evolved 
He donors(?). 



3.2. Properties of PLMXBs with Different Donors 

As fJT]mentions, the donors of LMXBs basically determine the orbital periods and mass-transfer 
rates which give the X-ray luminosity. Figure [2] shows the distributions of the orbital periods and the 
X-ray luminosities or the mass-accretion rate Mns- PLMXBs with different donors have different 
positions in Figure [2j 

? investigated the brightest Galactic PLMXBs which are plotted in Figure [2j They concluded 
that the majority of PLMXBs with X-ray luminosities below ~ 2 x 10^"^ erg s~^ have unevolved MS, 
while PLMXBs with higher X-ray luminosity predominantly harbor giant donors. In the panels of 
NS + MS in Figure O ~ 90% of PLMXBs with MS donors lie in the shallow region with a X-ray 
luminosity of ~ 10^^ erg s~^. The irradiation of LMXBs can drive mass transfer (??). In this work, 
we do not consider the effect of irradiation. Therefore, compared with the luminosities of known 
PLMXBs, we may underestimate the luminosity of PLMXBs with MS donors. In the panels of NS 
+ RG (Roche) in Figure [21 our sample covers the positions of Cyg X-2 and GX 13+1 which have 
long orbital periods. These PLMXBs can only be explained by NS + RG (Roche) systems in our 
simulations. ? gave good measurements for Cyg X-2 and have derived the donor's mass around 0.6 
Mq. ? suggested that the donor in Cyg X-2 has a mass of around 0.5 Mq with a non-degenerated 
helium core and is burning hydrogen in a shell. That is, the donor in Cyg X-2 is a sub-giant which 
has undergone the violent mass loss. This is consistent with ours. ? gave that the spectrum of 
donor in GX 13+1 clearly shows the features of the K5 III giant. Figure [3] gives the distributions 
of donors' luminosities and mass-transfer rates determine the X-ray luminosity. As the panels of 
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NS + RG (Roche) and NS + RG (Wind) in Figure [3] show, the donors' luminosities in PLMXBs 
with RG (Roche) donors are much lower than those in PLMXBs with RG (Wind), while the X-ray 
luminosities in the former are much higher than those in the later. Therefore, it is more difficult to 
observe the donors' luminosities in PLMXBs with RG (Roche) donors. We suggested that donors 
in Cyg X-2 and GX 13+1 are giants which fill up Roche lobes. The orbital periods and X-ray 
luminosities of 2 known symbiotic X-ray binaries (GX 1+4 and 4U 1700+24) are measured, which 
is plotted by triangles in Figure [2j NS + RG (Wind) systems can cover 4U 1700+24 very well, but 
our models can not explain GX 1+4. A detailed investigation of symbiotic X-ray binaries is being 
carried out (Lii et al. in preparation). 

In our simulations, the orbital periods of PLMXBs with WD or He donors are shorter than 
80 minutes, and they are ultra-compact X-ray binaries. About 10% — 34% of WDs in NS + WD 
systems are He WDs, and 66% — 90% are CO WDs. NS + ONeMg WD systems are negligible. 
As Figure [2] shows, PLMXBs with WD donors are agree with observations, while PLMXBs with 
He donors have X-ray luminosities higher than those of observed ultra-compact X-ray binaries. ? 
suggested that two ultra-compact X-ray binaries 4U 1626-67 and 4U 0614+09 have very evolved 
He donors, and their X-ray luminosities are 3.2 x 10^^ erg s~^ and 3.4 x 10^^ erg s~^. In our model, 
the mass transfer in PLMXBs with He donors is driven by the magnetic braking, which produces 
the X-ray luminosity of ~ 10^^ erg s~^. We may overestimate the work efficiency of the magnetic 
braking driving mass transfer. As panels of NS + WD and NS + He in Figure [3] show, the donors' 
luminosities in PLMXBs with He donors are much higher than those in PLMXBs with WD donors. 
This difference may be a way via which we can distinguish He donors from WD donors. 

Figure H] gives the distribution of the X-ray luminosities (mass-accretion rates) of NSs in 
PLMXBs with different donors. The X-ray luminosities between ~ 6 x 10^^ — 6 x 10^^ erg 
mainly come from PLMXBs with RG (Wind). The X-ray luminosities between ~ 10^^ — 10^^ 
erg s~^ originate from PLMXBs with MS, RG (Roche) and WD donors, in which PLMXBs with 
MS are dominated. The X-ray luminosities between ~ 10^^ — 10'^^ erg s^^ mainly originate from 
PLMXBs with He donors in which the mass transfer is driven by magnetic braking. ? suggested 
that observed X-ray luminosity function (> 2 x 10^'' erg s^^ ) of LMXBs can generally be explained 
by the accretion of matter onto a NS with magnetic stellar wind, which agrees with our results. 
However, the observed X-ray luminosity function ( up to ~ 2 x 10^'' erg s~^ ) of LMXBs can be 
explained by PLMXBs with MS and RG (Roche) donors. In BSE code, gravitational radiation is 
only efficient for binaries with orbital periods less than 3 hours. In order to explain a gap of CVs 
between 2 and 3 hr in the otherwise smooth period-mass distribution, BSE code does not apply 
magnetic braking when the primary is a fully convective MS star whose mass is lower than O.35M0. 
There is not the magnetic braking in NS + MS systems if MS's mass is larger than 1.2M0 because 
it has no convective envelope. In our work, less than 30% of PLMXBs with MS donors have a donor 
whose mass is between O.35M0 and 1.2M0. Therefore, the mass transfer in PLMXBs with X-ray 
luminosity function ( up to ~ 2 x 10^^ erg s~^ ) is mainly driven by stellar evolution or magnetic 
braking, but not gravitational radiation. 
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in the bottom panels of Figure HI we give the distributions of total X-ray luminosities which 
are different from the power-law X-ray luminosity function (showed by dashed line) observed by 
?. There are two main reasons: i)Our work does not include the transient X-ray luminosity. The 
transient systems in outburst reach high (close to Eddington) X-ray luminosities. ? suggested 
that the X-ray luminosities in outburst are between ~ 10^^ and 10^^ erg s~^. ii)We overestimate 
the X-ray luminosities of PLMXBs with He donors. Of course, we must obtain total population 
of LMXBs including transient and persistent systems with NSs and black hole before we compare 
X-ray luminosity function with observational one. In the future, we will try doing it. 



4. Conclusion 

We perform a detailed study of donors in PLMXBs, employing the population synthesis ap- 
proach to the evolution of binaries. We estimate that there are ~ 29000 — 110000 PLMXBs in 
the Galaxy, and their birthrates are ~ 3.4 — 7.2 xlO~'^ y^~^- PLMXBs with different donors have 
different formation channels. Our numerical simulation shows that more than 90% of PLMXBs 
have MS donors, and PLMXBs with RG (Wind) donors are neghgible. Most of NSs in PLMXBs 
with hydrogen-rich donors have undergone CCSN, while more than 90% of NSs in PLMXBs with 
He donors or WD donors have undergone AIC and EIC. 

In PLMXBs with MS donors, the orbital periods are between ~ 1 hour and 100 hours, and 
the mass transfer is driven by donor evolution or magnetic braking. Our population synthesis code 
shows that their X-ray luminosities mainly are around ~ 10^^ erg s~^. Similarly, in PLMXBs with 
RG donors via Roche lobe overflow(Roche), the mass transfer is driven by donor evolution, but 
orbital periods are between ~ 10 hours and 1000 hours. Their X-ray luminosities are around ~ 10^^ 
erg s-^. The 2 known LMXBs ( Cyg X-2 and GX 13+1) can belong to PLMXBs with RG (Roche). 
PLMXBs with RG (Wind) donors have the longest orbital periods and low X-ray luminosities. 
Their contributions to X-ray luminosities can be negligible. In PLMXBs with He donors, the 
orbital periods are shorter than 80 minutes, and the mass transfer is mainly driven by magnetic 
braking. Results of our numerical simulations predict that PLMXBs with X-ray luminosities around 
~ 10^^ erg s~^ mainly come from binaries with He donors, but their X-ray luminosities may be 
overestimated in our work. In PLMXBs with WD donors, the orbital periods are shorter than 1 
hour, and the mass transfer is mainly driven by gravitational radiation. According to results of 
our population synthesis code, their X-ray luminosities are between ~ 6 x 10^^ — 10'^^ erg s~^, and 
most of LMXBs with WD donors are transient. 

In this work we do not consider transient LMXBs with different donors. In further work, we 
will investigate donors in persistent and transient LMXBs, and discuss X-ray luminosity function 
of LMXBs. 
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Fig. 1. Number distribution of the masses of different donors in PLMXBs. The width of the 

bin is 0.005 Mq for PLMXBs with WD donors, and they are 0.5 Mq for others. 
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Fig. 2. Distributions of the orbital periods vs. the X-ray luminosities or the mass-accretion 

rate Mns in PLMXBs. The gradations of gray-scale correspond to the regions where the number 
density of systems is, respectively, within 1 - 1/2, 1/2 ~ 1/4, 1/4 ~ 1/8, 1/8 - of the maximum of 
giogP ^lagL ■ Number in every panel is normalized to 1. Circles represents the brightest PLMXBS 
from ?, and triangles represent the symbiotic X-ray binaries in ?. Solids represent the relation for 
hydrogen-rich disk of Eq.dS]), and dashed, dotted and dash-dotted lines represent the relations for 
He-rich, C-rich and 0-rich disk of Eq.dl]), respectively. 
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Fig. 3. — — Similar with Figure [2l but for distributions of donors' luminosity vs. the X-ray 
luminosities or the mass-accretion rate Mns in PLMXBs. 
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Fig. 4. Number distribution of the X-ray luminosities (mass-accretion rates) of NSs in PLMXBs 

with different donors. The width of the bin for logLx(erg/s) is 0.5. The dashed hnes in the total 
panels represent the power-law X-ray luminosity function for LMXBs in ?. 



